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ABSTRACT 

A new method arising from a gauge-theoretic approach to general relativity is applied 
to the formation of clusters in an expanding universe. The three cosmological models 
(f2o=l, riA=0), (rio=0.3, Qa=0-7), (rio=0.3, riA=0) are considered, which extends 
our previous application (Lasenby et al. 1999, Dabrowski et al. 1999). A simple initial 
velocity and density perturbation of finite extent is imposed at the epoch z = 1000 and 
we investigate the subsequent evolution of the density and velocity fields for clusters 
observed at redshifts z = 1, z = 2 and z = 3. Photon geodesies and redshifts are 
also calculated so that the Cosmic Microwave Background (CMB) anisotropics due 
to collapsing clusters can be estimated. We find that the central CMB temperature 
decrement is slightly stronger and extends to larger angular scales in the non-zero JIa 
case. This effect is strongly enhanced in the open case. Gravitational lensing effects 
are also considered and we apply our model to the reported microwave decrement 
observed towards the quasar pair PC 1643-1-4631 A&B. 

Key words: Gravitation - cosmology: theory - cosmology: gravitational lensing - 
cosmic microwave background - quasars: individual: PC1643-I-4631 A&B ~ galaxies: 
clustering 



1 SECONDARY GRAVITATIONAL EFFECT 

■ Rees and Sciama (1968) suggested that the presence of an 
evolving structure on the line of sight of a Cosmic Microwave 
. Background (CMB) photon could significantly affect its ob- 
' served temperature. This secondary gravitational effect is 
often described in terms of the potential well experienced by 
the CMB photon. For example, the potential well of a col- 
lapsing cluster becomes deeper over time so that the CMB 
photon has to climb out of a well deeper than that into which 
it fell, suffering a net loss of energy. As suggested in Rees and 
Sciama (1968) there is however a competing effect due to the 
extra time delay encountered by the photon. The overall ef- 
fect can therefore be of either sign. The photon accumulates 
a redshift along its geodesic, resulting in a net temperature 
perturbation. In the weak-field approximation we have (see 
Martinez-Gonzalez, Sanz & Silk 1990) 



that the photon traverses it, where <E>c is the gravitational 
potential of a rich Abell cluster. In this case, we have 

AT 



< 



y ^ (2) 

The potential of the cluster can be related to the velocity 
dispersion a through — GM/R = <I>c, where M and 
R axe the mass and radius of the cluster. Assuming a = 
gives 



AT 



< -2.2 X 10" 



(3) 



AT 



2 /• a$ , 



(1) 



where $ is the gravitational potential of the perturbation 
and t is the cosmic time. A rough estimate of an upper 
limit to this effect can be estimated by assuming that the 
potential well varies from $ = to "I> = $c during the time 
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In most cases, this anisotropy should therefore be small com- 
pared with other secondary anisotropics caused by the inter- 
action of CMB photons with non-linear structures, such as 
the thermal Sunyaev-Zel'dovich (SZ) effect. However, equa- 
tion (0) is only valid in the linear regime and a fully rel- 
ativistic treatment in the non-linear regime is necessary to 
estimate accurately the Rees-Sciama effect. 

Several methods have been used to model the formation 
of galaxy clusters in an expanding universe. The modelling of 
arbitrary perturbations in a fully general relativistic man- 
ner has not proved possible due to the complexity of the 
calculations involved. The most popular method has been 
the linearised approach (equation hi) , allowing the modelling 
of arbitrarily shaped inhomogeneities (Martinez-Gonzalez 
& Sanz 1990; Martmez-Gonzalez et al. 1990; Chodorowski 
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1991). This method has been particularly used for Great 
Attractor-hke structures located at redshifts as high as 
z = 5.9, where less strong non-linear effects are expected 
(Saez, Arnau & FuUana 1993, 1995; Arnau, FuUana & Saez 
1994; FuUana, Saez & Arnau 1994). 

In this paper we are interesting in modelhng rich galaxy 
clusters which are often ellipsoid or even quasi spherical (e.g. 
Coma cluster). This leads to spherical symmetry as a natu- 
ral approximation, which allows a general relativistic deriva- 
tion. The earliest attempts of modelling are of the Swtss 
Cheese type, where the collapsing Friedmann-Robertson- 
Walker (FRW) region is compensated for by a surrounding 
region of vacuum matched onto the expanding FRW uni- 
verse (Rees & Sciama 1968; Nottale 1982; Nottale 1984). 
In these models, analytic calculations can be performed in 
the fully non-linear regime and the observer can be located 
in the external, unperturbed universe. However, the matter 
distribution is unrealistic. 

More recent models have used the Tolman-Bondi so- 
lution for a pressureless, spherically symmetric matter dis- 
tribution to provide a more realistic density perturbation 
(Panek 1992; Quilis, Ibafiez & Saez 1995; Quilis and Saez 

1998) . These, however, are rarely compensated and comov- 
ing observers are not in a region modelled by a homogeneous 
FRW metric, which adds difficulty in estimating the temper- 
ature anisotropy. 

Finally, a new model arising from a gauge-theoretic ap- 
proach to gravity was proposed by Lasenby et al. (1999). 
This method employs a gauge, mainly characterised by a 
unique global time coordinate, which reduces the equations 
to an essentially Newtonian form. The model avoids the 
problem of streamline crossing, while keeping the density 
profile compensated and realistic. Furthermore, the pertur- 
bation is of finite size, allowing the observer to be placed in 
a non-perturbed region. 

In this paper we extend this new theoretical approach 
to the case of a non-zero cosmological constant A and to 
the case of open cosmologies. The non-zero A model is of 
particular interest since recent observations of high redshift 
Type la supernovae are inconsistent with A = fiat or open 
cosmologies, suggesting that our universe may be acceler- 
ating (Perlmutter et al. 1997, 1999). When the Type la 
supernovae results are combined with observations of the 
CMB, an approximately ffat FRW cosmological model is 
suggested with Qq ~ 0.3 and JIa « 0.7 (Efstathiou et al. 

1999) . fio is the mass density parameter at the present epoch 
and Ha = A/3Hq, where Ho is the Hubble constant at the 
present time. 



2 ASSUMPTIONS 

Three cosmological models are considered here: The Ein- 
stein de Sitter (EdS), Flat-A and Open models, correspond- 
ing to (^0 = 1, O,A=0), (f^o=0.3, Qa=0.7), and (^0=0. 3, 
r2A=0) respectively. Results will be quoted using the reduced 
Hubble parameter h{t), which is equal to H{t) expressed in 
units of ICQ km s~^Mpc~^. Furthermore, quantities evalu- 
ated at the present time t — to will be written simply with 
a zero subscript, e.g., h{to) — ho. Throughout this work we 
take ti, the time at which the initial perturbation is applied. 



to represent the epoch z = 1000. We employ natural units 
G = c = h = 1, unless stated otherwise. 

Throughout this paper, we assume that the cosmolog- 
ical ffuid as well as the evolving non-linear structure are 
pressureless. Quilis et al. (1995) included a hot gas compo- 
nent in their model and found that the pressure effects are 
in fact negligible and that the collapse is well approximated 
by the pressureless assumption. 

As mentioned above, the perturbation is assumed to be 
spherically symmetric which, combined with the pressure- 
less assumption, leads to the unrealistic situation where the 
structure collapses to form a singularity. As a result, the es- 
timates of CMB anisotropics given in this work should be 
regarded as upper limits. However, the initial perturbations 
are chosen so that the CMB photon traverses the cluster be- 
fore it becomes singular. Throughout this paper, the maxi- 
mum baryonic density experienced by the photon at the cen- 
tre of the non-linear structure is set to lO'^h^^^ protons m~"^. 
Furthermore the clusters considered in this work are located 
at high redshifts {z = 1, z = 2 and z — 3) and should there- 
fore still be in the process of formation. In these cases the 
large infalling motions predicted by the spherically symmet- 
ric pressureless assumptions may be less unrealistic. 

We further assume that, at all times, the baryon com- 
ponent contributes Wh^^^^ per cent of the total mass, the 
remainder being dark matter. The ft-dependence chosen here 
ensures that the total gravitational mass scales 
expected within the isothermal sphere approximation (e.g. 
Peebles 1993) . In this case the total density follows the power 
law form p cx and the baryon mass alone scales as hg^^^ . 



3 FLUID DYNAMICS 

This paper follows the theoretical approach of Lasenby et al. 
(1999) and Lasenby, Doran & Gull (1998) where exact and 
fully relativistic derivations are given using a new gauge- 
theoretic approach to gravity. The resulting field equations 
are remarkably straightforward and are here presented as 
if derived from classical first principles (see Lasenby et al. 
(1997) for rigorous calculations). 

The forming cluster and the external universe are both 
described by the same dynamical equations, the centre of 
the cluster being at the origin of the spherical coordinates 
(r, 6, (j)). Since we are only concerned with spherically sym- 
metric systems, we choose hereafter 6 — 7r/2 without loss of 
generality. The symbol t is used for the dynamical time. 

The two scalar quantities needed to describe the ffuid 
dynamics are the density p{t,r) and the radial velocity 

/ N dr ,., 
u{t,r)^— (4) 

along matter geodesies. The total gravitational mass M{t, r) 
contained within a sphere of radius r is defined by 

or 

M{t,r)= I 'iTvr^p{t,r')dr'. (6) 
Jo 

The mass of fluid flowing through a sphere of radius r per 
unit time is given by 
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— —4nr pu. 



(7) 



\ dt J r 

One can define the streamline derivative (i.e. co-moving with 
the fluid) as 

d 



_ = A A 

dl dt dr 



(8) 



We can see from ^ and (^) that the total gravitational 
mass enclosed is conserved along the fluid streamlines, 



^=0 

dl ' 

which forbids the possibility of streamline crossing. 



(9) 



3.1 Equation of Continuity 

By equating the total mass of fluid flowing out of a given 
volume per unit time to the decrease per unit time of fluid 
in the same volume, one obtains the standard equation of 
continuity 



where 



Hit,r) 



du{t, r) 
dr 



(10) 



(11) 



3.2 Euler's Equation 

Newton's second law applied on a fluid particle gives the 
usual Euler's equation 



(12) 



where F is the total force applied on the particle. F includes 
the radial gravitational force Fg — —M/r'^ as well as the 
repulsive force Fa proportional to r due to the cosmological 
constant A. To agree with the usual convention we set the 
constant of proportionality to A/3 so that Fa = (A/3)r. 
Finally equation (^2|) reads 

d. ... ^^3^ 



_ _du du _ M A 
dl ~~di^'^^^^7^^Y' 



3.3 Bernoulli's Equation 

Multiplying Euler's equation ( |l^ by u and integrating along 
the streamline dl, gives 

^"^-(v + r)^^(*'^)' (1^) 

where the constant of integration E{t, r) is constant along 
the fluid streamlines, i.e. 



dE 
"dT 



0. 



(15) 



In equation (p^), can be regarded as the kinetic energy 
of a fluid particle, while —M/r — (A/6)r^ can be thought of 
as its potential energy composed of the usual gravitational 
potential and a potential proportional to due to the pres- 
ence of a cosmological constant. E{t, r) can therefore be as- 
sociated with the total energy of a fluid shell. If S < the 
shell will eventually collapse at the origin while if i5 > 
the shell will remain in expansion. In other words, equation 



( p^ can be regarded as a generalised Friedman equation for 
which each fluid shell is allowed a different constant E. 



3.4 Newtonian Gauge 

The physical significance of the above set of equations is 
remarkably straightforward yet, as shown in Lasenby et al. 
(1998), they are fully consistent with general relativity. This 
is because the degrees of freedom offered by the gauge the- 
oretic approach have been carefully chosen so that the re- 
sulting equations look Newtonian. In particular the gauge 
fields have been fixed so that the coordinate r is a measur- 
able distance scale related to the strength of the tidal force, 
and the coordinate i is a measure of the proper time r for 
all observers comoving with the fluid, therefore dt = dl = 
dr, I being the afflne parameter of the streamlines (equa- 
tion Therefore this gauge choice enables us to recover a 
global "Newtonian" time on which all observers can agree. 
Lasenby et al. (1998) have named this gauge the Newtonian 
Gauge. The associated line element is given by 



ds' 



2E- 



dr^ + 



2u 



2E + 1 



-dtdr 



1 



2E + 1 



-dr 



-r^ (de'^ + sin^ 



(16) 



Similar analytic solutions have been found by Tolman (1934) 
and Bondi (1947), however their use of M(r) as the radial 
coordinate rather than r has the disadvantage of hiding the 
physical meaning of the field equation and complicates the 
choice of initial conditions. 

In homogeneous regions the fluid behaves as for stan- 
dard cosmologies and the density p is a function of time only, 
so we have 



47r Q 

M{t,r) = — r'p(t), 



(17) 



where the overbar denotes quantities in the external homo- 
geneous region. Here, equations (^, (|l^) and (^3|) lead to 



ft - -3^/''(^)' 
r 



and 
dH 



A 



(18) 
(19) 

(20) 



^ + ^ =-Y^"+3- 

It is clear from those standard results that we can identify 
H as the Hubble constant, and that the time parameter 
t of the Newtonian gauge agrees with the cosmic time in 
homogeneous regions of our model. We also have 



= H{tfm, 

for the EdS and 



en models and 



g-p(t) = H\t) - |A, 
for the Flat-A model. 



(21) 



(22) 



3.5 Streamline Equations 

In general, the solution to equation (^4|) can be expressed 
in terms of elliptic integrals. We choose here to integrate 
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numerically this equation using the following system of two 
first order ordinary differential equations composed of 
and Euler's equation: 

dr 

— — u 
(it 

du 

lit 



A 

— r, 
3 ' 



(23) 



where is the fluid particle radius on a given streamline, at 
an initial time ti. This system allows us to compute the po- 
sition r of the fluid particle as well as its velocity u at a later 
time t. We note that, in this manner, the change from pos- 
itive u to negative u is automatically evaluated in the case 
of closed streamlines. The fluid density p{t,r) and velocity 
gradient H(t, r) are calculated by differentiating numerically 
equations (g) and (pT[). 

However, equation (^) has analytical solutions in two 
particular cases, (i) For A = 0, which corresponds to the 
standard Friedman models. This case is treated in detail in 
Lasenby et al. (1999). (ii) For E = 0, where the solution is 
given by the following streamline equation 



t = ti 



+ 



^arccosh^ ^ 



/3A 
1 



\3M{r 



-r^ + 1 



arccosh ( „ . , , — ^rf -I- 1 



(24) 



V3A \SM{u 

In this case, the value of the velocity along the streamline is 
given by (yj) itself: 



2M{ri 



A 



(25) 



We note that here u is always positive since E — 0. Following 
Lasenby et al. (1999), the density is given by 

(26) 



where 



(f) 



Mr 2 y/Ar^ + 6M + 



6M 



(27) 



The velocity gradient is found by differentiating equa- 
tion (0), which leads to 



M 



dr 



+ A— + - ' 

6u u 



\drJt 



1 dM dE 
r dri dri 



(28) 



4 CLUSTER FORMATION MODEL 

The fluid equations given above control the evolution of the 
external homogeneous cosmology as well as the formation 
of the central cluster. The initial perturbation is applied at 
an early time ti that we choose here to represent the epoch 
z = 10^. The perturbation is of finite extent and we therefore 
need to fix initial conditions so that standard cosmology 
is satisfied outside a given radius FU. If ere we will use a 
family of simple four-parameter models based on polynomial 
perturbations in the density and velocity fields, as described 
in detail in Lasenby et al. (1999) and Dabrowski et al. (1999). 
The parameters are (i) the width of the perturbation Ri ; (ii) 



the velocity gradient at the origin H{ti, 0); (iii) the degree of 
the polynomial describing the perturbation, so that at r = 
Ri , the velocity and its first m derivatives match the external 
values; (iv) the external velocity gradient H{ti) equal to the 
Hubble constant at the time ti. As a typical value, we assume 
in this paper m = 3. 



4.1 Linearised Equations 

The perturbation responsible for the formation of clusters of 
galaxies is assumed to have grown from primordial fluctua- 
tions in the very early universe. At the epoch z = 10^, the 
fluid dynamics should still be well described by considering 
evolution in the linear regime. 

We introduce the linearised variables 



5u = u — rH(t) 
Sp = p- p(t) 



(29) 



5M 



A-KS 5p(t, s)ds , 



where the overbarred quantities denote background homo- 
geneous universe values. Differentiating equation (M) gives 



d&u 
"dT 



5M 



(30) 



as in the A = case (see Lasenby et al. 1999). Differentiating 
again and ignoring second order terms yields 



d /dSu 
di y'df 



- d(5it 



H^Su + ASu = 0, 



(31) 



where the background universe is assumed to be spatially 
flat {E = 0). By parameterising 5u in terms of streamline 
starting point, Vi and cosmic time, t, the streamline deriva- 
tives become derivatives with respect to t: 



Sil + 3HSu - {H^ - A)5u = , 



(32) 



where the overdot denotes derivatives with respect to t. In 
the limit of small t and A 



H : 



2 A 
3^+6^' 



and Su goes as 

5u^At~^^'' + Bt^^''+0{t''^^). 



(33) 



(34) 



Assuming that we are in an epoch where the decaying mode 
t"'*''"^ can be ignored but sufficiently early that the t'^^^ term 
is not yet significant, then 

(35) 



dt 



Alternatively, it can be said that at early enough times, 
2> A and so the A term may be ignored in equation ( ^2| ) . 
Thus using equation (^^ we have 

ISSU = -'-^, (36) 
2 

or expressed in terms of the perturbed and un-perturbed 
quantities 



u — rH 



2M r 87r_ 
3//r2 ^ 3H~S^' 



(37) 



Substituting 87rp/3 with values of equations ( pl[) and ( |22| ) 
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gives an extra condition on the initial data 

u(ti,r} — — ^ U 

for the EdS and Open models and 
2r /„ ,-,2 M 1 ' 



(38) 



(39) 



m \ r3 6 

for the Flat-A model. The initial density perturbation can 
now be derived from the initial velocity field as follows: Dif- 
ferentiating (pq) and (pfl) with respect to r gives 



p{ti,r) = — 



(3 + no)i?-^-(f^' 

r \ar y ti 



for the EdS and Open cosmologies, and 



3H 
8^ 

for the Flat-A case. 



2u 
r 



\drJt 



3H 



(40) 



(41) 




0.05 0.1 
radius (hp ^ Mpc) 

Figure 1. Initial relative velocity u{ti,r) — H(ti)r where ti rep- 
resents the epoch z = 10^. The solid, dashed and dotted lines are 
for the Open , Flat-A (Ca = 0-7) and EdS models respectively. 



4.2 Initial Conditions 

At the time ti, the initial conditions are defined completely 
by the 4 parameter velocity profile while the corresponding 
initial density profile is given by equation ( ^ ) or (^l[). The 
external universe H{ti) is fully defined by the desired cos- 
mology (i.e. f2o, f^A and ho). We have for the EdS a.nd Open 
models 



H{ti) = Ho{l + Zi)^/T+llozl, 
and for the Flat-A cosmology 



H{ti) = \Kl + z 



4- 



(42) 



(43) 



The two remaining parameters Ri and H{ti, 0) are fixed such 
that the resulting cluster is similar to those observed today, 
using the following characteristics: the redshift of the clus- 
ter Zc, its core radius Rc (where Rc is defined as the radius 
at which the cluster density falls to one-half its maximum 
value) and its central baryonic density pc. Because we are 
interested in modelling rich clusters, the core radius is fixed 
to Rc = 0.23/i(^^ Mpc and the central density is taken to be 
pc — W^Iiq^'^ p m~"^. The ft-dependences of Rc and pc are 
chosen so that the cluster's observed angular size and X-ray 
fiux are independent of ho (Jones & Forman 1984; Peebles 
1993) . Figs. |l| and |2| show the initial velocity and density pro- 
files in the Open , Flat-A and EdS models. In all three cases, 
the resulting cluster observed at Zc = 1 satisfies the same 



,^4 1,1/2 

10 «Q p m 



properties Rc — 0.23hg Mpc and pc 
We note that low-f^o models require larger and deeper per- 
turbation. This can be explained as follows: The equation 
equivalent to (^2|) but for the density perturbation evolu- 
tion is given by 



S = A-jvpS - 2HS, 



(44) 



where S — Sp/p. The first term of the right hand side gov- 
erns the growth of the perturbation while the second can 
be regarded as a damping term due to the drag from the 
expansion of the universe. In the low density universe cases 
(Flat-A and Open models), the growth term is smaller than 
in the critical density EdS model. Furthermore, H is larger 
in the Flat-A and Open cosmologies (see equations Ell and 




0.05 0.1 
radius (h^ ^ Mpc) 

Figure 2. As in Fig. ^ but for the initial density contrast 
p{ti,r)/p{ti). 



psj), enhancing the damping term. Therefore it is to be ex- 
pected that the initial perturbations in the case of the Flat- 
A and Open models have to be larger than in the EdS case 
in order to give rise to similar structures at a later time. 
A more detailed study shows that in the Flat-A model the 
growth suppression is less marked than for the Open model 
(see Peacock 1999). This explains why the initial perturba- 
tion of the Flat-A case is less pronounced than for the Open 
case, as seen in Figs. ^ and |^. 



4.3 Cluster Characteristics 

As discussed in Section ^, the pressureless and spherically 
symmetric assumptions inevitably lead to non linear struc- 
tures with large infall velocities, which would probably be 
realistic only for distant clusters still in the process of for- 
mation. Clusters observed at redshift z — 1, z = 2 and z — 3 
will be considered here and in the following sections. 

The central density profile obtain ed a t z — 1 when ap- 
plying the initial conditions of Section 4.2 is given in Fig. ^ 
up to the 1.5/i(^^ Mpc Abell radius, where the Open model is 
assumed. We found that the central density profiles obtained 
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0.5 1 1,5 

radius (hp ^ Mpc) 

Figure 3. The baryonic density of the cluster as a function of 
radius at the time corresponding to the redshift z = I. The dashed 
hne is the best-fitting /3-model up to r = 1.5/i^ Mpc. The Open 
model has been assumed here (see text). 



20 40 60 

radius (ti^ ^ Mpc) 

Figure 4. The velocity distribution relative to the Hubble flow 
u{t, r) — H{t)r as a function of radius at the time corresponding 
to the redshift z = 1. The solid, dashed and dotted lines are for 
the Open , Flat-A and EdS cosmologies respectively. 



in the context of the EdS and Flat-A cosmologies are almost 
indistinguishable from that of the Open case. Furthermore, 
varying the redshift at which the cluster is observed from 
2 = 1 to z = 3 does not affect the central density profile ei- 
ther. It is remarkable that by only fixing the desired central 
density pc and core radius Rc, the resulting density profile 
is almost identical in any cosmologies and at any redshift 
despite the significant difi'erence between the initial pertur- 
bations, as seen in Figs. |l] and |^. However, it was found in 
Dabrowski et al. (1999) that the shape of the density pro- 
file is sensitive to the other free parameter of our model m, 
which is here set to m = 3. 

As discussed in detail in Dabrowski et al. (1999), the 
obtained density profile is well fitted by /3-models 

Po 



p{r) = 



(45) 



usually used to fit observed profiles and rotation curves. It 
was also found that density distributions such as those de- 
rived from N-body codes (Navarro Frenk & White 1996, 
1997, hereafter NFW) give good agreement as regard the 
corresponding mass distribution of the cluster. However, the 
best least-square fit is obtained for a /3-model with /3 = 2/3. 
Indeed, as mentioned in Section ^ our model predicts density 
profiles p oc away from the origin which is more consis- 
tent with a, f3 — 2/3 model than with the NFW density 
profiles which behave like r'"*. Another difference between 
the density profiles obtained here and those of NFW is re- 
garding their behaviour at the centre of the cluster. In this 
work the central density keeps a finite value while it diverges 
to infinity in the case of NFW density profiles. 

The velocity distribution, as evolved at the time where 
the cluster is observed, is shown in Fig. ^ for z = 1 and 
for the three EdS , Flat-A and Open models. As expected 
from an initial perturbation of finite extent, the velocity 
profile matches exactly that of the external homogeneous 
universe after a finite radius. This effect can be seen clearly 
on Fig. ^ where the velocity relative to the Hubble fiow 
vanishes for r > 28hQ^ Mpc in the EdS cosmology and for 
r > 55hQ^ Mpc in the Flat-A and Open cases. It is not sur- 
prising that this radius is roughly a factor of 2 smaller in the 



Table 1. Cluster properties, z denotes the redshift at which the 
cluster is observed (see text). Rt is the turnaround radius quoted 
in hg^ Mpc and ( is the ratio of the central cluster density to that 
of the external Universe. Mi, Mi. 5, M2 and M4 are total gravita- 
tional masses contained within spheres of radii r = 1, 1.5, 2 and 



4/i„ Mpc respectively. All masses are quoted in 10 °h,, M 



10. 



Model 


z 


Rt 


c 


Ml 


Afl.5 


Ma 


M4 


EdS 


1 


6.26 


1122 


1.05 


1.80 


2.65 


6.70 


Flat-A 


1 


10.74 


3704 


1.00 


1.75 


2.50 


6.05 


Open 


1 


8.70 


3704 


1.00 


1.70 


2.50 


6.05 


EdS 


2 


3.22 


332.5 


1.05 


1.90 


2.90 


7.90 


Flat-A 


2 


6.01 


1097 


1.00 


1.85 


2.65 


6.75 


Open 


2 


4.94 


1097 


1.00 


1.85 


2.65 


6.75 


EdS 


3 


1.97 


140.3 


1.10 


2.10 


3.15 


9.50 


Flat-A 


3 


3.82 


463.0 


1.05 


1.90 


2.80 


7.50 


Open 


3 


3.23 


463.0 


1.05 


1.90 


2.80 


7.55 



EdS model since the initial velocity and density perturba- 
tions were smaller in the same proportion, as seen in Figs. ^ 
and H 

Other quantities of particular interest are the 
turnaround radius Rt , which is defined as the radius at which 
the fluid velocity changes from being radially inwards to ra- 
dially outwards, and the ratio ( of the cluster central den- 
sity to that of the external universe. These characteristics, 
together with total gravitational masses contained within 
spheres of given radii, are summarised in Table ^ for clus- 
ters at redshifts z = 1, z = 2 and z = 3 in the three EdS 
Flat-A and Open cosmologies. In all cases, the observed 

p m 



clusters have the same central density pc = 10 ilq 
and core radius Rc = 0.23/i(7^ Mpc. 



,4.1/2 „ „.-3 



5 PHOTON PATHS AND REDSHIFTS 

In order to estimate the effect of non linear structure forma- 
tion on the CMB we need to calculate the photon geodesies 
and redshifts in the context of the metric given in equation 
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Table 2. Estimates of the central temperature decrements AT/T 
in units of 10~^ for clusters at redshift 2 = 1,2 = 2 and 2 = 3. 
All clusters have the same central density pc = lO^hJ^^ p m~^ 
and core radius Rc = 0.23h^^ Mpc. Three cosmological models 
are considered. 



Model AT/T, 2 = 1 AT/T, 2 = 2 AT/T, 2 = 3 

EdS -0.64 -0.45 -0.32 

Flat-A -0.80 -0.64 -0.51 

Open -0.52 -0.43 -0.36 



Figure 5. Geometrical arrangement defining the photon geodesic 
integration. The grey discs denote the spherically symmetric col- 
lapsing cluster which is located at the origin of the coordinate 
system. The photon is emitted from a source co-moving with the 
Hubble flow. The photon path is defined by the photon (r, <j>) posi- 
tion as well as the direction x a^nd frequency ui of its 4-momentum 
P- 



(p^, as carried out in Lasenby et al. (1999). The photon mo- 
tion is parametri sed b y the Newtonian time parameter t dis- 
cussed in Section 3.4, The path is defined by four quantities, 
the photon radius r{t), position angle ^(t), 4-momentum di- 
rection x{i) ^''id frequency uj{t). As shown in Fig. H, our 



model allows us to start a photon path at a given source posi- 
tion co-moving with the Hubble flow; integrate the geodesic 
equations through the homogeneous universe and through 
the collapsing cluster; and finally collect the photon at the 
observer position, co-moving with the Hubble flow as well. 
The calculation is identical in any spherically symmetric and 
pressure free cosmological model so that the set of four first 
order differential equation given in Lasenby et al. (1999) is 
still valid here. 



6 EFFECT ON THE CMB 

As discussed in Section 0, the presence of a time dependent 
potential on the line of sight will affect our observation of the 
CMB temperature. For the model considered here, Lasenby 
et al. (1999) found 



AT r f dA 2 A . 2 



(46) 



where the integral is evaluated along the photon path. A is 
the fluid velocity relative to the Hubble flow 



A(t, 



rH{t), 



(47) 



ti is time at the surface of last scattering and t2 is the time 
at the observer. Numerical results are given in Table ^ for 
CMB photons that travel straight through the centre of col- 
lapsing clusters ai z — 1, z — 2 and z = 3. The EdS , and 
Open cosmological models considered here give roughly the 
same predictions for the central temperature decrement. We 
find however that the central anisotropy is slightly stronger 
in the Flat-A case. When looking at the redshift dependence, 
Table. ^ shows that, for a given cluster central density and 
core radius, the further away the cluster is, the smaller the 
magnitude of the Rees-Sciama effect is. At earlier epochs, 
collapsing structures are less decoupled from the Hubble ex- 
pansion and are therefore exposed to less strong non-linear 




-4,8 -4,7 -4,6 -4,5 

time (IIq ^ 10® years) 

Figure 7. Solid line: The integrand of equation (^), as a function 
of cosmic time t, along a photon path traversing a cluster located 
at 2 = 1 with an angle of view 9 = 0.5 degrees from the centre of 
the cluster. The time t=0 corresponds to the current epoch, and 
an Open model is assumed here. The central density of the cluster 
is Pc = W^hg^^ p m""^ and its core radius Rc = 0.23/Iq ^ Mpc. 
The corresponding temperature anisotropy AT/T is the area un- 
der the curve. The dashed and dotted lines are for (—dA/dr) and 
(— A/r) respectively. 



effects, therefore it is not surprising that the effect gets dim- 
mer. 

It is also useful to estimate the temperature anisotropy 
as a function of the projected angle 8 on the sky from the 
centre of the cluster. Profiles of AT/T are given in Fig. ^for 
the three EdS , Flat-A and Open models and for clusters at 
2 = 1,2 = 2 and 2 = 3. As noticed earlier, the magnitude of 
the effect gets smaller as the redshift of the cluster increases. 
Fig. ^ reveals that the angular distribution of the anisotropy 
is strongly dependent upon the cosmological model. In par- 
ticular for the case where the cluster is located at 2 = 1, the 
Open model estimate predicts a large and extended temper- 
ature increment on a scale of about one degree. 

The fact that the temperature anisotropy can take posi- 
tive as well as negative values is not surprising when looking 
at equation (0). As long as 8 is not too small {—dA/dr) is 
negative, while (—A/r) is always positive, as seen in Fig.^. 
Typically, there are a negative and a positive contributions 
to the temperature anisotropy AT/T which are weighted by 
cos^ X and sin'^ x respectively. For larger 8, sin^ x is bigger 
and the positive exceeds the negative contribution. The in- 
tegrand of equation ( ^ij ) is plotted in Fig. |^ along the path 
of a photon passing near the centre of a cluster located at 
2 = 1. The positive and negative contributions to the re- 
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0,5 1 1.5 2 10 0.2 0,4 0.6 O.B 1 10 0.2 0.4 0.6 



d (degrees) 9 (degrees) 9 (degrees) 

Figure 6. Estimates of the angular distribution of the CMB anisotropy due to the presence of a collapsing cluster on the line of sight. 
From left to right, the three figures are for clusters at redshift z = 1, z = 2 and z = 3 respectively. All clusters have the same central 
density pc = lO'^h^^^ p m~^ and core radius Rc = 0.23h^^ Mpc. The anisotropy is in units of AT/T as a function of the projected angle 
6 on the sky from the centre of the cluster. The solid, dashed and dotted lines are for the Open , Flat-A and EdS cosmological models 
respectively. Note that the angular scale is different for each figure. 



suiting AT/T are separated out. In this case, the overall 
perturbation is AT/T ~ +0.17 x 10"^ 

In the case of the recently favoured Flat-A cosmological 
model, it is interesting to note from Fig.^ that the temper- 
ature decrement is systematically stronger and extended to 
larger angles than in the other two models. 



7 GRAVITATIONAL LENSING AND TIME 
DELAY 

In this section we apply our model to the problem of 
the microwave decrement reported towards the quasar pair 
PC 1643+4631 A&B (Jones et al. 1997). This decrement 
is observed at 15 GHz and lies between a pair of quasars 
at redshifts z = 3.79 and 3.83. The two quasars are sepa- 
rated by 198 arcsec. If the microwave decrement is assumed 
to be a thermal SZ effect caused by an intervening cluster, 
the minimum central temperature decrement is estimated to 
be Ar/r=— 2.1 X 10~^. However deep optical and infrared 
imaging carried with the WHT and the UKIRT (Saunders 
et al. 1997, Haynes et al. 1999) as well as ROSAT X-ray 
observations (Kneissl, Sunyaev & White 1998) suggest that 
any intervening cluster producing the temperature decre- 
ment should be located at a redshift as high as z = 3. How- 
ever, the presence of a cluster is supported by the remark- 
able similarity between the quasar spectra, suggesting grav- 
itational lensing. Saunders et al. (1997) suggest that their 
~ 1 per cent redshift difference can be explained in terms 
of quasar intrinsic evolution over a delay between the two 
lightpaths of ~ 10'^ years. It was found in Dabrowski et al. 
(1999) that the time delays obtained in the EdS cosmology 
were too short. Here we propose to carry out a similar inves- 
tigation but in the Flat-A and Open cosmological models. 



Table 3. Gravitational lensing time delays in the frame of the 
emitting source located at z = 3.8. The lens is placed at a reshift 
2 = 1, 2 = 2 or 2 = 3. Rc is for the core radius of the cluster 
quoted in Mpc and At is the time delay in years. In all cases 
the cluster has a central density of pc = lO'^h^^^ p at the 

time photons crosses it and the lens system is characterised by an 
Einstein ring of 100 arcsec. 





2 = 


1 


2 = 


2 


2 = 


3 


Model 


Rc 


At 


Rc 


At 


Rc 


At 


EdS 


0.45 


145 


0.75 


205 


2.68 


25 


Flat-A 


0.40 


250 


0.60 


550 


1.36 


380 


Open 


0.44 


195 


0.70 


425 


1.89 


140 



In this section, we assume Ho = 50 km s~^Mpc~^; the 
maximum baryonic density along the photon path is still 
fixed to 10** p m~'^; the total gravitational mass over bary- 
onic mass ratio is 10. Following Dabrowski et al. (1999), we 
simplify the situation by considering a source directly behind 
the centre of the cluster. This source of photons is co-moving 
with the Hubble flow and is located at a redshift z — 3.8. The 
cluster, located at z = 1, z = 2 or z = 3 is the gravitational 
lens and, in order to be consistent with Jones et al. (1997) 
we require the Einstein ring of the lens system to be 100 arc- 
sec. Such a constraint allows us to fix the core radius Rc of 
the cluster, as given in Table |^. The estimated time delays 
At are between a photon from the source travelling straight 
through the centre of the cluster and one which follows a 
lensed path, appearing at the Einstein ring radius. Results 
are given in Table ^ as a function of the cosmological model 
assumed and the redshift of the cluster. Since the source is 
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at fixed redshift z = 3.8 the lensing effect has to be stronger 
for more distant clusters in order to keep an Einstein ring 
of 100 arcsec. This is why the core radii are getting larger 
with redshift, as seen of Table |^. Consequently, the time de- 
lay is increasing as well. In the EdS model, this is probably 
leading to rather un-realistic situations such as for z = 3 
where the core radius of the cluster is as high as 2.68 Mpc 
(369 arcsec). In this case, the light path passing through the 
Einstein ring at 100 arcsec is also travelling through the core 
of the cluster. This is why the estimated time delay of 25 
years is so small. We find that in the Open and particularly 
in the Flat- A cosmological models, the predicted time de- 
lays are larger than for the EdS case and may be enough to 
explain the observations related to PC1643-(-4631 A&B, if 
they are indeed a gravitational lens system. 



8 CONCLUSION 

We extend a general relativistic model for the formation of 
non-linear cosmic structures (Lasenby et al. 1999; Dabrowski 
et al. 1999) to the case of open cosmologies and non-zero 
cosmological constant. The equations of fluid dynamics are 
presented as if derived from first principles. This is possi- 
ble when working in the context of the Newtonian gauge 
(Lasenby et al. 1998). The fluid streamlines are derived an- 
alytically in the case of a flat universe with non-zero cos- 
mological constant. However, in non- homogeneous regions 
the streamlines have to be computed numerically. The ini- 
tial conditions for the fluid at the early epoch correspond- 
ing to z = 1000 are very simple: a flnite perturbation on the 
fluid velocity is imposed, and the corresponding density per- 
turbation is inferred assuming a linearised model and that 
the perturbation arose from primordial fluctuations. As ex- 
pected, in order to give rise to similar structures at a later 
time, we find that the initial perturbation has to be larger in 
the case of a low density universe than in a critical density 
universe. 

Photon geodesies can be integrated numerically (see 
Lasenby et al. 1999) and we study the gravitational effect of 
the collapsing galaxy clusters on CMB photons. Clusters at 
redshifts of z = 1, z = 2 and z — 3 are considered. We find 
that the central temperature decrement AT/T is systemat- 
ically slightly stronger and extended to larger angles in the 
non-zero cosmological constant case. It was also found in 
the open cosmology model that the CMB anisotropy could 
reach relatively large positive values. 

Finally, we apply our model to the microwave decre- 
ment reported towards the quasar pair PC 1643-1-4631 A&B 
(Jones et al. 1997). Here we are particularly interested in 
estimating the time delay in the source frame between the 
two quasar images. We find a time delay of ~ 550 years in 
the case of a fiat universe with non-zero cosmological con- 
stant (Jlo = 0.3, Qa = 0.7). If the quasar images are indeed 
lensed images from a single source, this time delay may be 
sufficiently large to explain the 1 per cent difference in the 
quasar redshifts. 
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